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ABSTRACT 

We present an in-depth study of the High Mass X-ray Binary (HMXB) GX 301-2 during its pre-periastron flare using data from the 
XMM-Newton satellite. The energy spectrum shows a power law continuum absorbed by a large equivalent hydrogen column on the 
order of 10^"* cm"^ and a prominent Fe Ka fluorescent emission line. Besides the Fe Ka line, evidence for Fe K/3, Ni Ka, Ni K/3, 
S Ka, Ar Ka, Ca Ka, and Cr Ko- fluorescent lines is found. The observed line strengths are consistent with fluorescence in a cold 
absorber. This is the first time that Cr Ka is seen in emission in the X-ray spectrum of a HMXB. In addition to the modulation by 
the strong pulse period of ~685 sec the source is highly variable and shows different states of activity. We perform time-resolved as 
well as pulse-to-pulse resolved spectroscopy to investigate differences between these states of activity. We find that fluorescent line 
fluxes are strongly variable and generally follow the overall flux. The A'h value is variable by a factor of 2, but not correlated to 
continuum normalization. We find an interval of low flux in the light curve in which the pulsations cease almost completely, without 
any indication of an increasing absorption column. We investigate this dip in detail and argue that it is most likely that during the 
dip the accretion ceased and the afterglow of the fluorescent iron accounted for the main portion of the X-ray flux. A similar dip was 
found earlier in RXTE data, and we compare our findings to these results. 

Key words, stars: neutron (GX 301-2) - X-rays: binaries - Accretion, accretion disks 
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J — 1 . Introduction The exact type of the companion and its distance was loi ig dis - 

pO , „ , ,. ^ , TT- , T,- puted. but the most recent measurements bv lKaper et al.l (l2006h 

^ ^"JC-o^^^J^.n/^ coordinates of the High Mass X-ray Binary .^^^^^^ ^ ^.^^^^^^ ^^^^ .^^^^^ ^ ^ gl ^ 

O (HMXB) GX 301-2 was Al^etected in 1969 during a balloon -^^^ ^-^^^ ^ stellar-radius of 62 R/and a mass of 43 ± 10 m'. 

^ . experiment ILewin et a L 097^. During further balloon flights ^^^p^ these val ues throughout the paper, as X-ray accretion 

^ . l McClintocketal.1 im3 found that the source was persistent, ^^^^j^ bv lL^^fe Kostka (2008.) also seem to indicate similar 



\l ' Since then GX 301-2 has been observed by most X-ray instru- y^J^gg 
^ , ments and has stirred a lot of interest in the community as it 

k> shows a strongly variable flux and a very large absorption col- The regular and bright X-ray flares of GX 301 -2 were found 

/S umn. The source shows an erratical l y chan ging pulse period of to occur ~ 1 .4 d before the periastron passage, i.e., before the neu- 

^ ~ 685 sec, detected by IWhite et all (119761) and in the last few tron star reaches the densest pait of the stellar wind. Standard, 

■ ■ years monitored by the Fermi/GBM Pulsar Project ( Tinger et al.L symmetric wind accretion models fail to explain pre-periastron 

I2OO9E clearly identifying the compact object as a neutron star. fare. Pravdo & Ghosh ( 2001 ) suggested that a circumstellar disk 

The system has an orbital period of P^a = 41.506 + 0.003 d is formed around Wray 977, from which the neutron star ac- 

at flie reference time ro(MJD) = 43906.06 + 0.1 1 and the or- cretes additional matter when it passes through. If the disk is not 

bital period is cha i iging with Porb = (-3.7 + 0.5) x 10 Ss' aligned with the orbital plane of the neutron star, a pre-periastron 

(iDoroshenko et all 1^^). These values were determined by ac- Aare could be observed. More recent calculations (Leahyl |2002|; 

curat e pulse timing an d are in agreement with earlier results Leahy & Kostka, 2008) show that it is also possible that the neu- 

from lKohetall (Il997h who used the regular X-ray outburst of tron star is trailed by an accretion stream, which it overtakes 

GX 301-2. The eccentricity of the orbit was determined to be shordy before periastron. This model self-consistendy explains 

e - 0.47. the observed dependence of flux and column density on the or- 

The neutron star accretes matter from the optical com panion , bital phase. It is the regular pre-periasti-on flare which distin- 

which was identified as the B-type giant Wray 977 (,VidalllT97l . guishes GX 301-2 from most other HMXB. Du ring the flare 

the X-ray flux increases by up to a factor of 25 dPravdo et all 

' see also |http : //www . bat se . msf c . nasa . gov/gbm/science/pul"s ^i^/g;:iRih tfaachitA /fa>St(ft1irijai M HT i. reaching luminosities of up to 
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Fig. 1. S ketch of th e orbit of GX 301-2, using the parameters 
given bv iKohetalJ dl997h . The green line shows the average 
RXTE- ASM light curve as function of orbital phase. The phase 
of the XMM-Newton observation discussed here is highlighted 
in red. 



1 Crab in the 15-50keV band and some lOOmCrab in the softer 
2-10 keV band. X-rays with energies below ~5 keV are less vari- 
able over the orbit as they are strongly affec ted by scattering 
and a bsorption due to the high column density (iLeahv & KostkaL 
I2OO8 I). 

Figure □ shows the orbit of GX 301-2 and the RXTE- ASM 
light curve fo lded onto an orbi t al pe riod of 41.5 d using the 
ephemeris by iDoroshenko et al.l (I201Q) . Clearly visible is the 
pre-periastron flare, peaking around phase tp - 0. 91. The plot 
assum es a size of 62 for the optical companion (iKaper et al.L 
l2006h . 

The broad-band X-ray spectrum of GX 301-2 can be de- 
scribed by a power law with a high-energy cutoff around 20 ke V 
and a cyclotron re sonant scattering feature (CRSF) around 
37keV ("End o et al.L l2002h . lKrevkenl3ohm et al.1 (|2004") showed 
that the spectral parameters, e.g., the energy of the CRSF, vary 
with pulse phase, and therefore concluded that the X-rays are 
most likely produced in an accretion mound close to the neutron 
star surface. In the softer X-rays a very strong Fe Ka fluorescent 
emission line is visible together with a photo-electric absorp- 



tion with an equiva l ent hy drogen column A^h of up to 10 

(iLeahvet all 119881119891) . Using Chandra data . IWatanabe et al.l 



(1200 3h detected a fully resolved Compton shoulder at the low- 
energy flank of the Fe Ko' line. These authors could model the 
shoulder by simulating the Compton scattering of photons in a 
dense, cold gas-stream. 

The slow pulse period of GX 301-2 is strongly variable 
over long timescales ('Pravd o & GhoshL 1200 it lEvangelista et al.L 
L2OIO). It was first detected with P ~ 700 sec, but changed in 
the last 40 years to values as low as 675 sec and is today around 
P ~ 685 sec. Although the period is variable, GX 301-2 usu- 
ally shows clear p ulsations throughout its orbit. Recently, how- 
ever, iGogii^ et al.l (1201 II) . discovered a dip in the RXTE X-ray 
light curve during which GX 301-2 ceased to show pulsations. 
The data showed that exactly one pulse was missing and that the 
pulsed fraction was gradually reduced before and gradually in- 



creased after the dip. Interestingly, the spectrum softened signif- 
icantly during that dip, contradicting the common picture that a 
very dense blob, ab sorbin g all soft X -rays, was passing through 
the line of sight. iGogii^ etsU (|201 ih speculate that due to the 
cessation of accretion the black body radiation from the polar 
cap became visible, effectively softening the spectrum and re- 
ducing the pulsed fraction, as the polar cap is larger and un- 
dergoes stronger gravitational light bending than the accretion 
column. Similar di ps have already been obse rved in the light 
curves of Vela X- 1 dKrevkenbohm et ani2008h . but never before 
in GX 301-2. 

The strong absorption during the bright pre-periastron flare 
makes GX 301-2 an ideal object to study the X-ray producing 
and absorbing regions around a neutron star, especially focusing 
on the many fluorescent lines. XMM-Newton is the ideal instru- 
ment for this observation, with its large effective area and good 
spectral resolution. We obtained a ~47 ksec XMM-Newton ob- 
servation of which the results are presented here. 

In the next Section we will describe the data used. Section 3 
will cover the timing analysis and in Sect. 4 the different spectral 
analyses will be discussed. In Sect. 5 the results are discussed 
and a conclusion and an outlook are given. 



2. Data & Observations 

We present data taken with the XMM-Newton observa- 
tory on July 12th, 2009 (MJD 55024.103-55024.643, ObsID 
0555200401) at orbital phase -0.91 of GX 301-2, i.e., during 
the pre-periastron flare. The exposure was 46.6 ksec, covering 
only a small part of the flare (see Fig. [TJ. The X-ray flux was 
strongly variable during the observation, with values between 
3.67 X 10"'*' erg sec"' cm"^ and 2.47 x 10"^ erg sec"' cm"^ in the 
2-lOkeV energy band. Due to the expected high flux, the EPIC- 
MOS detectors were turned off to provide more telemetry for 
the EPIC-pn camera, which was operated in timing mode. In 
this mode the CCD is continuously read out, so that all spatial 
information in the read out direction is lost. In the other direc- 
tion only the central 64 columns are used. The timing mode al- 
lows the measurement of countrates up to 400 cts sec"' without 
detectable pile-up (Struder et al., 2001). As GX 301-2 shows 
very strong absorption no usable information can be gained from 
the Reflection Grating Spectrometers (RGS), which are mainly 
sensitive below 2.5 keV. We used the standard Science Analysis 
System (SAS) 10.0.0 to extract light curves and spectra. We 
tested for pile-up using the tool epatplot, but found no severe 
pile-up. The source data were extracted from EPIC-pn columns 
35^0 and background data were obtained from columns 10-17 
and 53-60. As high photon fluxes change the charge transfer in- 
efficiency (CTI) of the CCD detector, we also extracted spectra 
with the tool epfast, which takes these changes into account. 
We found, however, no significant differences to the standard 
pipeline. The data were barycentric and binary corrected, using 
the ephemeris provided bv .Doroshenko et al.. C2010.) . 



3. Timing Analysis 

3.1. Light curve 

The top panel of Fig. |2]shows the 20 sec resolution light curve of 
GX 301-2. The light curve shows clear flaring behavior, with a 
smaller flare in the beginning and a large flare starting ~30 ksec 
after the start of the observation. Clearly visible are also the sub- 
peaks induced by the pulse period of ~685 sec. We performed 
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Fig. 2. a) Lightcurve of GX 301-2 with 20 sec resolution in 
the 0.5-lO.OkeV energy band. Vertical lines indicate the time 
regimes I-IV used in the later analysis, b) Hardness ratio be- 
tween the continuum energy band C (8.5-lO.OkeV) and the ab- 
sorbed energy band A (4.0-5.0keV). c) Hardness ratio between 
C and the Fe Ka line band F (6.3-6.5 keV). 



a search for the pulse period via epoch folding and found a pe- 
riod of P = 685.0 + 0.1 sec, consistent with the monitoring data 
from the Fermi/GBM Pulsar Project. In search for changes of the 
pulse period during the observation, we performed epoch fold- 
ing and phase connecting on parts of the light curve as short as 
15 ksec. We did not find a significant change of the pulse period, 
consistent with Fermi/GBM data which shows that the pulse pe- 
riod was only slowly changing during that time. 

A close-up view of the light curve reveals that albeit the pul- 
sations are clearly visible, strong pulse-to-pulse variations are 
also present. The individual pulse shapes range from simple two- 
peaked profiles to more complex shapes, with multiple minor 
peaks superimposed on the general shape. An example of two 
different profiles is shown in Fig. fusing light curves with 1 sec 
and 20 sec resolution. It becomes clear that the average pulse 
profile, also shown in Fig. [3] is highly complex. Additionally 
brightness changes independent of the pulses are visible in the 
light curve. These changes can be on the same order with respect 
to luminosity and duration as the variability due to the pulsation. 

Stronger and sudden changes in luminosity are also seen in 
the light curve and different luminosity levels are marked as re- 
gions I-IV in Fig. |2] These changes are also accompanied by 
changes in the pulse profile and the pulsed fraction. To investi- 
gate this in more detail, we folded each regime I-IV separately 
on the pulse period of P = 685.0 sec and studied the normalized 
distribution of the countrate in each phase bin. As the overall 
countrate is different for every region we also normalized each 
light curve to its maximum countrate. The mean of the distribu- 
tion in each phase bin gives the value of the pulse profile in that 
bin. Figure |4] shows the results for all four regimes in a color- 
coded map. Region III shows the clearest pulse profile, where 
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Fig. 3. Two examples of the variability of the pulse shape, left 
an example early in the observation, right one close to the end 
of the observation. Shown in black is the 1 sec light curve, in red 
a 20 sec light curve. The blue profile is the average pulse profile 
over the complete observation. Note that the y-axes are not the 
same scale. 
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Fig. 4. Color-coded maps of the pulse phase distribution of nor- 
malized countrates in regions I-IV as indicated in Fig. |2l See 
text for details. Color represents the number of measurements 
for each countrate in the respective phase bin. Each phase bin 
was normalized to itself, so that they have all comparable color- 
codes. Colors are ranging from violet, corresponding to zero, to 
red, corresponding to one. The superimposed black hne shows 
the mean value, i.e., the pulse profile. 



the mean countrate as well as the whole distribution is changing 
with phase. This behavior is different in region I. Here the pulse 
profile is also visible in the mean values but the values measured 
most often are only weakly phase dependent. The two peaks of 
the profile result from strong outliers and flares at the respective 
phases. In region II a similar effect is seen, but the overall pulsa- 
tions are even weaker. Region IV has a very short duration and 
therefore provides barely enough statistics for an interpretation. 
The behavior seems similar to the one seen in region III, but the 
pulse profile is not that well pronounced. 

Region II seems to be the region with the weakest pulsations 
and shows a very different behavior than region III. To investi- 
gate these differences further, we took a closer look at the 1 sec 
light curve during region II, as shown in Fig. |5] In the beginning 
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Fig. 5. Small part of the light curve between 0.3-12.0keV, show- 
ing a vanishing of the pulsations. In black the light curve 
with 5 sec resolution is shown, in red with 40 sec resolution. 
Superimposed in blue is the average pulse profile for the whole 
observation. 

the pulse is visible but it vanishes almost completely during the 
decline of the luminosity for ~3 rotations. Around t - 21000 sec 
a sudden increase of the luminosity is accompanied by weak 
pulses but these disappear again after ~A rotations of the neu- 
tron star, also during a dip in the brightness. Looking at the hard- 
ness ratio during that time (Fig. |2] lower two panels) allows us 
to deduce the spectral shape. The hardness ratio was calculated 
as HR(//, S)-{H- S)/(H + S), with H indicating the count rate 
in the harder and S in the softer energy band. When comparing 
the count rate in the 8.5-lOkeV continuum band C to the count 
rate in the strongly absorbed energy band A between 4-5 keV, 
no significant variability of the hardness HR(C, A) can be mea- 
sured during region II. Neither is a change compared to region 
I visible, indicating that the reduced flux is not due to increased 
absorption, which would primarily absorb the soft X-rays and 
would thereby make the spectrum overall harder. Contrary to 
that we see that the hardness HR(C, F) between the count rate 
in the continuum band and in the iron line energy band F be- 
tween 6.3-6.5 keV changes dramatically in region II. The strong 
drop of the hardness indicates that the flux of the iron line is not 
reduced as strongly as the continuum flux. 

3.2. Pulse profiles and pulsed fraction 

In the energy range of XMM-Newton, the average pulse profile 
is only weakly energy dependent (Fig. |6]l. At higher energies 
the smaller, secondary peak becomes relatively stronger com- 
pared to the first one, but still stays clearly weaker. Only at en- 
ergies above 20keV the secondary peak gets almost as bright 
as the primary one (Kreykenbohm et al., 2004; La Barber aet al.L 
I2OOI. This behavior is also confirmed in recent Suzaku data 



dSuchv etal.L l2011h. 
Ilr 



llndo et al.l (|2002|) and lTashiro et alj (1199 lb have shown fliat 
the pulsed fraction in the iron line energy band is significantly 
smaller than in other neighboring energy bands. Analysis of 
XMM-Newton data confirms this behavior, indicating that the 
matter producing the fluorescent iron line is distributed homo- 
geneously around the X-ray producing region. The pulsed frac- 
tion / was calculated using the formula 



/ 



max(P) - min(P) 
max(P) + min(P) 



(1) 



with P denoting the pulse profile. The pulsed fraction is plotted 
in Fig. |7] for the four different regions defined above. Common 
to all four regions is the fact that the pulsations in the iron line 




Fig. 6. Pulse profiles in different energy bands: a) 0.3^.0 keV, b) 
4.0-6.2keV, c) 6.2-6.6keV, and d) 6.6-15.0keV. Each profile is 
shown twice for clarity. 
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Fig. 7. Pulsed fraction for different energy bands and for the 4 
different regions I-IV. The dashed vertical line marks 6.4 keV, 
the energy of the Fe Ka line, the dotted one 7.1 keV, the energy 
of theFe Ky61ine. 



band are notably reduced. In region IV the dip is least prominent, 
probably due to the overall low pulsations and the short duration 
of this region. 

Figure I2] also clearly shows that the pulsed fraction is vari- 
able with time, dropping significantly during region II, as already 
seen in Fig.|4] Some residual pulsed fraction is still measurable 
outside of the iron line band, but as the iron line band domi- 
nates the whole spectrum during that interval, almost no pulsa- 
tions are visible in the broad band light curve. The lowest energy 
bin, between 0.5^.0keV is dominated by the strong absorption, 
smearing out the pulsations and reducing the signal-to-noise ra- 
tio (S/N). With increasing energy the pulsed fraction increases 
above 6.4 keV, with exception of the energy bin covering the 
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Fe K/S line at 7.1 keV. The Fe KjS line is strong enough to con- 
tribute significantly to the flux in that bin, and the line is, as the 
Fe Ka hne, only very weakly pulsed. 

4. Spectral analysis 

4.1. Phase averaged spectrum 

When looking at the averaged X-ray spectrum of GX 301-2 
(Fig.|8}j) one sees a large iron Ka fluorescent emission line and 
an equally large absorption column, which reduces the countrate 
to effectively zero below 2 keV. The spectrum can be reasonably 
well modeled with a partially covered power law, with additional 
Gaussian lines at the Fe Ka, Fe KJ3, Ni Ka, Ca Ka, Ar Ka, and 
S Ka fluorescent line energies. We used a model of the form: 



<x 



^(cA^H.i +(1 - c)A^H,2)x power law + Gaussian lines j (2) 



Where A^h.i and A^h,2 are the two different parts of the partial 
coverer and < c < 1 is the covering fraction. We applied the 
tbnew absorption model (Wilms et al., 2011, in prep. Jj for both 
elements of the partial coverer. For the large columns measured 
in GX 301-2 the plasma becomes optically thick and some flux 
is scattered out of the line of sight. This effect is taken into ac- 
count with the cabs model, described by 



cabs(£') = exp(-A^H o-j(E)) 



(3) 



where crj(E) is the Thomson cross section and A^h is the equiv- 
alent hydrogen column. Both components of the partial cov- 
erer were multiplied by the cabs model using the same col- 
umn densities, A^h,i and A^h,2, respectively. The Gaussian com- 
ponents used to model the fluorescent lines were only photo- 
electrically absorbed using a galactic Nf^ = 1.7 x lO^^ cm 
(iKalberla et al.L l2005|) . i.e., they were modeled as they appear 
outside the absorberQ This description was chosen because the 
spectrum shows clear Ar and S Ka lines at 2.3 and 2.9 keV, re- 
spectively, which can not be as strong as observed if they have 
to pass through an absorption column of ~ 10^"^ cm"^, as is 
necessary to describe the continuum. For a discussion on self- 
absorption and intrinsic strength of the lines, see Sect. 15.21 This 
spectral model, including the same emission lines, can also de- 
scribe the broad- band spectrum of G X 301-2, as seen in re- 
cent Suzaku data (ISuchv et al.L 1201 ih . Spectral parameters can 
be found in Table [T]and a plot of the spectrum with the separate 
components and the residuals in Fig.|8] In Table[T]A^H,i gives the 
absorption column which is multiplied by the covering fraction 
(CF), while Nh,2 is multiplied by (1 - CF). The normalization 
and the photon index of the power law continuum are given in 
rows "norm" and "photon index F", respectively. 

The Fe K a hne shows strong r esiduals when fitted with a sin- 
gle Gaussian. IWatanabe et al.l (I2003.) have shown that the Fe Ka 
line shows an extended Compton shoulder (CS) at its lower en- 
ergy side, but the resolution of EPIC-pn does not allow one to 
identify this feature uniquely. Bianchi et al. (2002) have shown 
that the Compton shoulder can be described with a box-shaped 
model. We included such a box in our model, which extends the 
Gaussian shaped Fe Ka line to lower energies and starts at the 
centroid energy of the line. This model gives a very good descrip- 
tion of the line, as evident in Fig.|8] In Table[T]the row "CS Ef 



gives the low energy end of the box, while "CS A" gives its area. 
The ratio of the flux between the Fe Ka line and the Compton 
shoulder is in good agreement with the ratio expected for pho- 
tons transmitt ed through a dense absorption column on the order 
of 10^3 cm-2 (lMati[200l . 

The best fit obtained with this model was 769 with 
376 degrees of freedom (dof) Otred" 2.04). Still snake-like 
residuals between 5.0-5.5 keV were visible, see panel b) of 
Fig. [8l The Kg line energy of neutral chromium i s at 5.4 1 keV 
dKaastra & Mewel 11993^ iKortright & Thompson! l2001h . and 
adding a Gaussian emission line in this energy range let to an 
improvement of the fit to;if^ =713 with 374 d.o.f. The center of 
the line was found at 5.425^QQ3QkeV in agreement with theory. 
Still residuals remain around 5.6 keV in absorption, see panel c) 
of Fig. [8] These residuals are in the energy range of the plateau 
between the iron lines and their escape peaks in the silicon of the 
detector The escape peaks and the plateau are clearly visible in 
the model of the Gaussian components of the Fe Ka and Fe K/3 
line in Fig. [8] due to the large flux of these lines. The residuals 
could be very well described with a broad Gaussian absorption 
feature, improving the fit further to;^'^ = 488 with 371 d.o.f. and 
called "Abs. Feat." in Table[T] Using only this absorption feature, 
and leaving out the Cr Ka line, to model the residuals does not 
adequately describe the spectrum, and gives ax^ of only 562 for 
373 d.o.f., see Fig. [g/. 

Finally strong residuals remained around 8.3 k eV, very close 
to the energy of the Nickel K^ line at 8.26 keV (iPalmeri et al.L 
l2008al) . A Gaussian shaped line centered at that energy was used 
to model the residuals and improved the fit clearly. This last ad- 
dition resulted in our best fit model, with;^'- - 448 for 369 d.o.f. 
(X^^^ - 1.214), of which the residuals are shown in Fig.[8f! 

To our knowledge, this is the first time that clear evidence for 
a Cr Kff as well as a Ni KJ3 emission line was found in the soft 
X-ray spectrum of a HMXB. Cr has been observed in absorption 
in oth er sources before, e.g., in GRO J1655-40 ( Miller et al.L 
1 20081) . These authors found, however, no evidence for fluores- 
cence lines in the spectrum of GRO J 1655 -40 and Cr as well as 
all other Fe-like elements were found to be highly ionized. 

Even though the spectrum was only used above 2keV, as 
all bins below are clearly background dominated, the continuum 
in the energy range between 2-3 ke V can not be accurately de- 
scribed by only the partial covering model. A constant DC level 
is needed additionally, to lift the continuum in the lowest energy 
bins. Without this DC component the best fit provides only ax^- 
value of 737 with 370 d.o.f. Values of the DC level are given 
in Table [T] in the row labeled "DC level". The DC level might 
originate from a low energy response of the large Fe Ka line be- 
low the escape peak, which is not fully taken into account in the 
RMF A similar efi'ect is known to be present in XIS data from 
Suzaku dMatsumoto et al.Ll2006l ). 



4.2. Time resolved spectra 

We have seen in the variability of the light curves, pulsed frac- 
tions, and hardness ratios that the overall spectrum will only 
provide a very inadequate description of the physics of the neu- 
tron star, as many different accretion regimes are mixed together 
Sacrificing S/N we therefore decided to split the spectrum ac- 
cording to the four regimes I-IV defined in Sect.[3]and indicated 
in Fig.|2] These states show the strongest differences in the light 
curve and hardness ratio and are therefore expected to show dif- 



^ see also http : / /pulsar . sternwarte . u ni-erlangen.de/wilms/i6^ fe^B^QiMbs/ 

3 DL value obtained with nh (HEASOFT v6.1Q), provided by AU spectra were rebinned to a S/N of 12 up to 3.5 keV, of 36 



HEASARC, at the coordinates of GX 301-2 



up to 6keV, of 60 up to 8.5 keV and of 48 above. Additionally, 
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Energy [keV] 

Fig. 8. a): Phase and time averaged spectrum in the 2-lOkeV 
energy range of EPIC-pn. The components of the spectral model 
are shown in different colors, see text for a detailed description. 
b) Residuals with neither an absorption feature nor an Cr Ka line 
around 5.4 keV. c) Residuals with a Cr Ka line, but without an 
absorption feature, d) Residuals with an absorption feature, but 
without a Cr Ka line, e) Residuals without a Ni K/3 line./j Best 
fit residuals. 

at least 4 channels were combined to account for the real energy 
resolution of the instrument. This binning resulted in different 
amounts of degrees of freedom for each spectrum. To each of 
these spectra we fitted the same model as to the overall spectrum 
and obtained very good values ~1 for all four spectra. 

The best fit parameters are listed in Table [T] It can be seen 
that the absorption is highest during part II in which the coun- 
trate drops very low and the pulsations vanish. During that time 
the absorption is increased by ~2Q% compared to part I. The un- 
derlying power law continuum softens at the same time, result- 
ing in comparable hardness ratios. The unabsorbed flux is also 
highest during part II (see row iTj™^^ j^^y in Table [T]), so that the 
observed drop in count rate can be attributed solely to the very 
high absorption column. 



The width, cr, of the Fe Ka was in all parts around 40 eV, 
larger than the energy resolution of EPIC-pn. If this width 
is due to Doppler shifts from turbulence in the fluorescent 
medium, speeds in excess of 1000 km s"' would be necessary. 
In GX 301-2, observational evidence and models suggest that a 
wind speed o f only 300km s~' i s reached at the p osition of the 
neuti-on star (iKaper et all l2006t Il7eahv & KostkaL 120081) . which 
can account for significantly lesser broadening of only ~10eV. 
During accretion the wind will become turbulent, however, and 
thereby possibly locally increasing the wind speed to values 
high enough to explain the broadening dMauche et al.L l2008l) . 
Turbulent flows in the accretion stream are, however, theoreti- 
cally not yet well understood. A more likely source of apparent 
line broadening comes from the fact that the Fe Ka is not a sin- 
gle line, but a superposition of different Ka lines of differently 
strongly ionized iron (Kallman et al., 2004). This effect can ac- 
count for a width of 20-30 eV of the observed line, as the single 
lines cannot be resolved with the spectral resolution of EPIC- 
pn. Together with broadening in the wind, the observed width of 
Fe Ka seems realistic. The Fe Kj8 line is even broader with val- 
ues around 65 e V. As the Fe KB com plex is not as broad as the 
Fe Ka complex dKaflman et al. . 120041) . this broadening might be 
mainly due to confusion with the strong Fe K-edge at 7.1 keV, 
very close to the line. EPIC-pn's resolution does not allow to 
completely separate these two features, but attempts to fit the 
spectra with a narrow Fe K/3 line lead to unacceptable;(f^gjValues. 

The centroid energy of the Fe Ka line is changing with time 
and can be as high as 6.45 keV in part III, where the brightest 
flare occurred. This energy corresponds to Fe xix or higher ion- 
ized iron, which does not emit a Fe K/3 line, as all M-shell elec- 
trons are already stripped from the ion. Still the spectrum of part 
III shows clear evidence for a F e KJ3 line at 7.127ke V, corre- 
sponding roughly to Fe xiii to xiv (iKallman et al.Ll200 4^). For the 
Ni Ka and Ni KJ3 lines a similar discrepancy emerges regard- 
ing the best fit values: the best fit Ni Kfi c entroid energy is very 
close to the one of neutral nickel ( Palmeri et a l UIOOSa) . while 
the Ni Ka centroid energy requires at least Ni X or higher. The 
uncertainties of the Ni K/3 energy, however, allow also for ion- 
ization states of this order 

Regarding the flux ratios between the Fe Ka and Fe KJ3 we 
obtain strongly varying values between Fe K/3/Fe Ka = 0.186 - 
0.283. The Ni K/J/Ni Ka ratio seems to be less variable with 
values between 0.18 - 0.23, and the individual uncertainties 
are much larger so that they are all consistent with each other 
(see Table [T] rows "Ratio Fe" and "Ratio Ni"). These values 
are clearly above the theoretical ones for neutral elements from 
iKaastra & Mewd (l993), which give Fe K/3/Fe Ka = 0.125 and 
Ni KyS/Ni KJ3 - 0.105. Measurements in solid state metals give 
higher nu mbers of Fe K / 3/Fe K a = 0.132 and Ni K/3/Ni Ka 
= 0.133 dHan & Demiii l2009l) . The Ni K/3/Ni Ka ratio is 
marginally consistent with all these values within its uncertain- 
ties, but the measured iron ratio is clearly higher Other measure- 
ments and calculations for neutral iron give fig ures of the same 
order with Fe K/3/Fe Ka ~ 0.13 (see, e.g.. , Palmeri et al.Ll2003L 
and references therein). 

It is conceivable that the flourescence lines originate from 
the dense absorbing material. If this is the case, self-absorption 
is not negligible (llnouelll985h . For an column density of A^h ~ 
140 X 10^^ cm"^ the optical depth at the energy of the Fe Ka line 
is only =s 60% of the depth at the energy of the Fe KJ3 line. If 
the whole absorption column would act on the lines, the intrin- 
sic ratios are drastically reduced to values between Fe KJ3fFe Ka 
- 0.114 - 0.173, much closer to the theoretical values. Such a 
strong effect is, however, unlikely, as the fluorescence lines at 
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Table 1. Fit parameters for the best fit model for the overall spectrum and the spectra of the four states I-IV, respectively. For details 
on the model see text. Uncertainties are given at the 90% limit. The two lowest values are obtained when using a gainshift kernel, 
as discussed in the text. 



Parameter 



all 



I 



n 



m 



IV 



Coy. Frac. 
A'h.i [10^" atoms cm"-] 
Nn,2 [10^- atoms cm"-] 
Norm [keV"' s"' cm"-] 
Photo Index F 
Fe Kq-A [ph s"' cm"-] 
FeKaE [keV] 
FeKo-cr [keV] 
CS El [keV] 
CS A[phs"icm"2] 
FeKjSA [phs"icm"-] 
FeKjSE [keV] 
FeKj8o-[keV] 
Ni KaA [phs"' cm"^] 
NiKorE [keV] 
NiKc!'0-[keV] 
CaKcAEph s"'cm"2] 
CaKaE [keV] 
Ar Ka A [ph s"' cm"^] 
ArKo-E [keV] 
S Kq-A [phs"' cm"-] 
S KaE [keV] 
Abs. Feat. A [phs"' cm"^] 
Abs. Feat. E [keV] 
Abs. Feat, cr [keV] 
CrKffA [phs"'cm"2] 
CrKttE [keV] 
NiK/?A [phs"'cm"2] 
NiK/?E [keV] 
DC Level [keV"' s"' cm"-] 
Ratio Fe 
Ratio Ni 
Fe Ions 

T2_iokcV [erg sec ' cm -] 

mabs 



141+3 
53.7 ± 1.6 

1 1 c+0.00 
^•^-^-0.06 

0.90 + 0.02 

(i.04!»:«^] X 10"^ 

g4227+o.6oi6 

0.0426!»oo,8 
6.205 ± 0.009 

(2.0 ± 0.3) X 10"3 
(2.32 ± 0.05) X 10"3 
7.095 ± 0.002 
0.0675!»:™^ 
(0.57!»:«1) X 10"3 
7.491 ±0.005 
036+" °'" 
(0.73 ± 0.06) X 10"^ 

T 71+0.02 

■'• ' ^-0.01 

(3.2 ± 0.3) X 10"5 

9 Q7T+0.012 
^•^'-^-0.002 

(3.3 ± 0.4) X 10"^ 
9 ^i<;+0.008 

(-0.52:°:?^) X 10"' 

5.61 ±0.03 

0.26 ± 0.03 

(1.5 ± 0.2) X 10"'' 

5 430+""'^ 
-'•^-^ -0.000 

(1.2 ± 0.3) X lO""* 
8.26!SS 



943+" '" 
147 ±6 

59 ±4 

1 iq+0.18 
'^■'^'^-0.15 



0.93 ± 0.04 

(o-82!jJ:!!^)xio" 



6.413 : 
0.040 : 



0.003 
0.004 



6 199+0''" 

(2.4+;;;^) X 10"3 

(1.76 ± 0.08) X 10"3 
7.087 ± 0.005 
0.069 ± 0.006 

(4.4+«5) X 10"* 

7 4C4+0.005 

< 0.04 

(0.47 ± 0.09) X 10"* 

T 7 1 +0.03 
' ^-0.00 

(1.9 ± 0.5) X 10"5 

9 QfiQ+O.Oie 
■^•^"^-0.014 

(2.1 ±0.5) X 10"^ 

2.35 ± 0.02 

(-3.3:] ^) X 10"* 
c C7+0.06 

(1.2 ± 0.5) X 10"* 
c 45+0.04 
-^■^-^-0.02 

(0.8 ± 0.5) X 10"* 

Q 00+0.06 
O-^°-0.08 



(0.48 ± 0.04) X 10"* (3.9 ± 0.6) x 10" 



r,^!'-^l^ [erg sec"' cm"^] 



fl 99^+0007 
"•^^-^-0.006 

0.20 ± 0.05 



448.02/369 (1.214) 
1.043 X 10"' 
1.617 X 10"** 



214+0014 

"■"^^^-0.013 

0.18 ±0.11 
XII-XIII 

113.69/122 (0.932) 
8.863 X 10"'o 
1.624 X 10"** 



9796+00"'^ 

"■^'^"-0.0018 

175!^ 

47 ±2 
9 9+0.3 
^•^-0.1 
1.10 ±0.04 

(I.10t°l)xl0"^ 

6.407 ± 0.002 

0.040 ± 0.003 

(2.2+0*) 10-3 

(2.04 ± 0.07) X 10"3 
7 070+0 003 

'■"'"-0.004 

0.063 ± 0.005 

(0.66!0-05) X 10"' 

7 470+0.004 
'■^'-^-0.000 

053+00" 
"■"-'-^-0.012 

(0.82 ± 0.09) X 10"* 

3 716+o oot 
-'•' ^"-0.011 

(3.1 ±0.5)x 10"5 

9 Q5C+0.015 

Z..>'JJ_QQQQ 

(3.0 ± 0.5) X lO""' 

2.315 ±0.011 

(-2.3!? *) X 10"* 

5 67+0 03 
•'■"'-0.02 

12+0'' 
(0.9!0]) X 10"* 
5.42 ± 0.03 
(1.2 ± 0.4) X 10"* 

g 99+0.03 
"■^^-0.04 

(3.4 ± 0.6) X lO""' 

'J-ioO-o.oov 
0.18 ±0.06 
VI~XI 

136.13/120 (1.134) 
7.127 X 10"'o 
2.151 X 10"** 



89+0 0^ 

"■"'^-0.03 

114±5 
53+3 

1 02+009 
0.78 ± 0.03 
(l.27!0-05)xlO"2 

6.453 ± 0.003 

032+0 003 

6.249ii 
(3.3!-)xl0"3 

(3.68!0-li] X 10"3 
7 l2669+"^ oooo9 

'•^^""'-0.00016 

0.052 ± 0.005 
(0.69 ± 0.07) X 10"3 
7 515+0-015 
<0.03 

(1.03 ± 0.18) X 10"* 

T 7QC + O.OI5 

(0.47 ± 0.09) X 10"* 

3.000 ± 0.015 

(0.61 ±0.10) X 10"* 

2.348 ± 0.012 

(-1.5±0.5)x 10"3 

5 70+0 00 
-OQS 
45+0 09 
"•^■^-0.12 

(l.9!0^)x 10"* 

5.44 ± 0.03 

(1.9 ± 0.8) X 10"* 

o 97+0.05 
°-^'-0.04 



+ 1.2 ■ 



0.968! 
130: 
42:f 

1 05+008 



(0.78!0-?9) X 10"2 

6.417+0-008 

029+^-^''^ 

(l.3t2)xlO"3 
(l.53!0-|5) X 10"3 
7.063 ± 0.009 
0.058!0:0|5 
(0.42!0-|5) X 10"3 

7.50 ± 0"03 
06+0 05 

"•""-0.06 

(0.8 ± 0.2) x 10"* 
3 69+0 02 

■^•"'^-0.03 

(0.48^0 1^) X 10"* 
2.96 ± 0.02 
(2.9+[-^) X 10"^ 
2.36 ± 0.06 
(-0.8!03) X 10"3 
5 56+o 'o 

-DOS 

32+0-'^ 
"•-'^-0.12 

(1.9 ± 0.7) X 10"* 

s 44+0.02 
-^■^^-0.03 

(0.8!0;^)xlO"* 
9 9+0.3 



(0.74 ± 0. 1 1) X 10"* (0.4+0 5) X 10" 



n 90 7+0.020 

"■^'^■'-0.016 

0.23 ±0.14 
>XIV 

129.76/131 (0.991) 
1.852 X 10"' 
1.797 X 10"** 



20+0 04 
"■"^"-0.02 
0.2 ±0.2 

IX-XI 

91.72/90 (1.019) 
8.048 X 10"'o 
1.350 X 10"** 



Gainshift 
^Vd.o.f(y?„,) 



1.00396 ±0.00017 
468.08/375 (1.248) 



1 0023+"-™ 



1 0010+"™3 

^■""^"-0.0002 



129.48/128 (1.012) 172.64/126 (1.370) 



1.0082 ±0.0003 
133.93/137 (0.978) 



1 0019+"™^ 
101.12/96 (1.053) 



softer energies from elements like S and Ar would vanish com- 
pletely. The line ratio also depends strongly on the viewi ng angle 
and the geometry of the absorber dMolendi et al.L l2003l and ref- 
erences therein). Simulations have shown that for certain view- 
ing angles around 45° of a dense slab absorber in transmission 
the measured ratios are reduced, while softer fluorescence lines, 
especially the Cr Ka line, are still clearly detected (Barragan et 
al., 2011, in prep.). A detailed investigation will be presented in 
a forthcoming publication. For further discussion of the absorp- 
tion of the fluorescence lines see also Sect. 15.21 

The line energies as well as the flux ratio is expected to 
change with the ionization state of iron (P almeri et al., 2003). 
Allowing for a mixture of different ioniza tions during each par t 
and using the energies and yields from Mendoza et al.l (l2004l) . 
we found that ions between Fe vi to Fe xiv can explain the ob- 
served energies and ratios (see Table[T] row "Fe Ions"). Here we 
neglected self-absorption, as we can not confine its effect on the 
ratios. Additionally the Fe Kfi line is superimposed on the iron 
K-edge and the Gaussian model component can be artificially 



increased due to a coiTelation with the edge. The real ionization 
states are therefore likely to be lower, but changes are still clearly 
present. We find an indication that the ionization state is corre- 
lated with the luminosity of the source but that overall the ion- 
ization state stays around intermediate values. It is furthermore 
possible that the energy calibration of XMM-Newton is slightly 
off and that a small shift in energy depending on the count-rate 
is present, leading to an apparent shift of the observed line en- 
ergies. The most likely effect causing such a behavior is either 
a change in the charge transfer inefficiency (CTI) of the CCD 
or a gainshift due to X-ray loading. The CTI is decreasing with 
higher ph oton flux when impu rities in the CCD are filled up with 
electrons (IStruderetal.Lll997i) . while X -ray loading causes more 
electrons than expected. If these changes are not taken into ac- 
count correctly, they lead to a shift of measured energies depend- 
ing on the source flux. The expected uncertainties of the calibra- 
tion of these effects in EPIC-pn are below +0.5% (Guainazzi, 
2011, priv. comm.). 
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To investigate the possibility of such a shift we froze the line 
energies of all measured elements to the values of their respec- 
tive neutral line energy, with values taken from Palmeri et aD 
(2008b) for S Ar, and Ca, [Kaastra & Mewe (1993) for Cr, 
fKallman et alj (|2004|) for Fe, and Palmeri et al.i (i2008a) for Ni. 
We then applied a gainshift kernel, shifting the model counts 
spectrum energy grid according to £'new(c) = iSoidCc) /gainshift. 
The obtained best fits were only marginally worse than with free 
energies of the lines and the shifts were on the order of 1.001- 
1.008 (0.1-0.8%), see bottom Unes of Tableffl These values are 
only marginally higher than expected from calibration uncertain- 
ties but also include unknown source intrinsic energy shifts. The 
continuum parameters of the spectral model did not change sig- 
nificantly due to the gainshfit, but the fluxes of the fluorescence 
lines were slightly altered to values marginally consistent with 
the ones shown in Table[T] Especially the Fe K/S line flux was de- 
creased, possibly due to a better seperation from the iron K-edge. 
This effect influenced also the Fe KfifFe Ka ratio by decreasing 
it slightly to values <0.2 at all times, which is closer to the values 
expected for neutral iron. No significant change with luminosity 
could be observed. 

Using the gainshift kernel we can explain the measured line 
energy shifts stemming purely from detector effects while the 
physical condition of the plasma in GX 301-2 does not change 
significantly during our observation. 

4.3. Pulse to pulse spectroscopy 

In an approach to understand the time variability of the spec- 
tral parameters, especially the iron line complex, in more detail, 
we extracted spectra for individual pulses, leading to 69 spectra 
with 685.0 sec exposure time each. The resulting model param- 
eters are presented in Figs. |9]and[T0] The data are grouped by 
color into 9 parts labeled A-I to facilitate orientation. We ap- 
plied a reduced version of the time averaged spectral model and 
included only the Fe Ka line with its Compton shoulder, the 
Fe KJ3, and the Ni Ka line. All other fluorescent lines could not 
be identified in these short exposure spectra. Additionally we 
had to fix the photon index F - 0.90 and the secondary absorp- 
tion column A^H,2 = 5.37 x lO^"* cm"^ because neither could be 
constrained well with the data. When leaving the centroid energy 
of the emission lines free, we found a strong correlation between 
these energies and the average count rate, as already indicated in 
the previous Section. Th is correlation is clearly stronger than ex- 
pected from ionization dKallman et al.L l2004l) and is therefore 
likely a detector effect. As described in the previous Section, 
these calibration uncertainties can be modeled with a gainshift 
kernel. We also checked the spectra extracted using epfast, but 
they showed the same systematic effects. 

With fixed energies for the fluorescent lines and for the iron 
K-edge the use of a gainshift kernel leads to a very good de- 
scription of the individual spectra. As expected the gainshift is 
directly correlated to the incident countrate, see Fig.fTOfe. 

The flux of the Fe Ka line and its Compton shoulder is 
mostly, but not always, correlated with the count rate in the 0.5- 
lOkeV energy band, as seen in Fig.|9] Especially in the second, 
main flare (part G) the iron line flux starts to rise with the light 
curve but breaks down around the maximum of the flare and 
stays low for the remainder of the observation. A better corre- 
lation is observed with the normalization of the power law con- 
tinuum (Fig. [Tol l. The Fe KJ3 and the Ni Ka lines show a very 
similar behavior and are strongly coiTelated with the Fe Ka line 
flux, with a Pearson's correlation coefficient of 0.92 and 0.89, re- 
spectively. As discussed in the previous Section, the Fe Ka and 




Time [sec] 

Fig. 9. Evolution of the spectral parameters of individual pulses. 
a) Light curve between 0.5-10 keV , b) Fe Ka flux (including 
the Compton shoulder), c) Fe KJ3 flux, d) Ni Ka flux, e) ratio 
between the Compton shoulder and the Gaussian shaped Fe Ka 
line, and f) ratio between the Fe Ka and the Fe Kfi flux. The 
fluxes of the lines are given in photons s"' cm"^ x 10"^. 



Fe Kj3 lines required a finite width cr, with an average value of 
40 eV and 71 eV, respectively. The widths did not show a signif- 
icant correlation with time or flux. The width of the Ni Ka was 
frozen to 10 '*keV. 

The absorption column is also strongly changing over the 
course of the observation, varying by more than 50% (see 
Fig.fTOb). It is increasing during the first flare in part C and stays 
high until the onset of the second, main flare (parts F and G), 
where a strong drop is registered. After the flare it is increas- 
ing again gradually. The covering fraction is very high through- 
out the observation, with an average of ~92%, reducing the sec- 
ondary absorption column to almost negligible values. It is only 
during the bright flare in part G where the covering fraction sud- 
denly drops as low as 74%, around the time when the Fe Ka line 
flux dips. 

A typical spectrum with one pulse exposure time is shown in 
Fig HH The Fe Ka and Fe KJ3 lines as well as the iron K-edge, 
are clearly seen. The Ni Ka line is also visible on close inspec- 
tion, but below 3 keV no significant excess over the background 
can be measured anymore due to the huge absorption column. 
The Ca, Ar, S, and Cr Ka lines and the Ni K^ line are not signif- 
icantly detected. 
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10000 20000 30000 40000 

Time [sec] 

Fig. 10. Time series of the spectral parameters of the spectra of 
individual pulses, a) Light curve between 0.5-lOkeV, b) con- 
tinuum norm in photons keV"' s ' cm c) main absorption col- 
umn in 10^^ cm~^, d) covering fraction, e) gainshift slope, and/) 



-2 ^of the best fit. 



Assuming that the absorbing matter is the same matter re- 
sponsible for the fluorescent iron emission and that it is dis- 
tributed spherically symmetric around the X-ray source, the 
equivalent width of t he Fe Kg lin e is proportional to the over- 
all absorbing column dlnou Furthermore, self-shielding 
should be negligible and the matter is assumed to be cold and 
neutral. We can see this behavior during most of the XMM- 
Newton -observation, even at very large A^h values, see Fig. [T2l 
The average Pearson's correlation coefficient over all parts ex- 
cept part D is 0.91. During part D, corresponding to the off-state 
seen in the light curve (Figs.|2]and|5]l, this correlation is not vis- 
ible. Instead indication for an anti-correlation with a correlation 
coefficient of -0.67 exists. During that state the equivalent width 
of the Fe Ka line reaches its highest values, whereas the A^h 
stay s in a regime around 1.8x10^"* cm"^, lower than the maximal 
measured value. The overall A^h value is quite high, as expected 
from the analysis of region II from Sect. 14.21 (corresponding to 
parts D, E, and F) which showed the highest overall absorption 
value. But the off-state is happening during part D only, a part 
which does not stand out for its increased absorption column 
when compared to parts E and F. This invariance of the A^h value 
clearly shows that the off-state is not due to a increased absorp- 
tion, which was already indicated by the softening during that 
state (see Fig.|2|. 




Energy [keV] 

Fig. 11. Typical spectrum with 685.0 sec exposure time, i.e., one 
pulse. The upper panel shows the data and the best fit spectrum, 
the lower panel shows the respective residuals. For details about 
the spectral model see text. 
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Fig. 12. Equivalent width of the Fe Ka line as function of the 
main absorption column A^h.i- Colors are the same as in Fig.fTOl 
and are indicating the time. 



5. Discussion & Conclusion 

We have presented a detailed timing and spectral analysis of the 
HMXB GX 301-2 during its pre-periastron flare using the EPIC- 
pn camera of XMM-Newton. The 2-lOkeV spectrum can be de- 
scribed with a partially covered, strongly absorbed power law 
model enhanced by several Gaussian emission lines from S, Ar, 
Ca, Fe, and Ni. We found evidence that the Ka line from Cr is 
also present as well as the K/? line from Ni. To our knowledge 
this is the first detection of these lines in a HMXB. In order to ac- 
count for the different states of activity visible in the light curve 
we separated the data according to these states, labeling them 
I-IV. During the sub-flares (parts I and III) regular pulsations 
with P = 685.0 sec were clearly detected and a double peaked 
pulse profile is visible. During the off-state II these pulsations 
vanished almost completely and the countrate dropped to very 
low values of ~10% of the peak count rate. XMM-Newton still 
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detected the source clearly. By separating the data into individ- 
ual spectra for every rotation of the neutron star, we were able 
to perform a detailed spectral analysis of the off-state. We found 
that although the absorbing column is very large, it is not dras- 
tically increased during that state, compared to values measured 
shortly before and after it. Instead a softening of the spectra is 
observed in which the Fe Ka line is even more important than at 
other times. 

We have shown that the absorption column is larger than 
10^"^ cm"^ at all times during the observation. GX 301-2 
is not the only source showing absorption columns of this 
magnitude. INTEGRAL, for example, discovered more than 
10 highly absorbed HMXB, which might for m a new class 
of objects (iKuulkersl 120051; IWalter et al.[ l2006h. Among theni 
are sources like IGR J16318-4848 (iMatt & GuainazziL 
I2OOI .B arragan et al., 2009) or IGR J17252-3616 



dManousakis & Walteii 1201 II) which show very prominent 
Fe Ka and Fe Kfi lines, as well as Ni Ka lines. The flux of 
these sources is, however, usually lower than that of GX 301-2, 
decreasing the S/N drastically below 4keV and thereby not 
allowing for studies of softer fluorescence lines. Results from 
the iron line energies and Fe Ka/K/3 ratios also indicate only 
mildly ionized iron. GX 301-2 could therefore be a brighter 
example of these sources and help to understand this class of 
HMXB. 

5.1. Fluorescence region size 

To estimate the distance between the X-ray producing region and 
the fluorescent medium, we extracted light curves with 1 sec res- 
olution for the iron line energy band between 6.3-6.5 keV and 
the harder continuum band between 7.3-8.5 keV. The countrate 
between 6.3-6.5 keV is completely dominated by the Fe Ka line 
and any variability of it must be due to the reaction of the flu- 
orescent line to changes in the continuum flux. We compared 
both light curves using Fourier transformation as well as cross- 
correlation techniques and could not find a significant time delay 
between both curves in the range of 2-5000 sec. This is in agree- 
ment with the size of the fluoresce nce region being smaller than 
10^" cm, i.e., 0.3 It-sec, as given bv lEndo et alJ (12002) . The count 
rate in the narrow energy bands was not high enough to achieve 
a time resolution sufficient to see a del ay of the fluorescence 
lines on that order. For IGR J17252-3616 Manousakis & Walter' 
(1201 1 ) also found that the emission line region is compact and 
not larger than the accretion radius. However, during the off- 
state in part D in Fig. |9] the iron line is still clearly visible, 
although the luminosity of the continuum has significantly de- 
creased. This leads to the largest equivalent widths for the Fe Ka 
line measured in the whole observation (Fig. [T2l l and seems to 
indicate that the fluorescence line is partly reacting rather slowly 
on changes in the continuum. A possible explanation is that the 
weaker part of the partial coverer is located further out, and 
thereby reacting much slower to any changes in the X-ray contin- 
uum flux. The X-ray flux was low for ~4500 sec and we can as- 
sume this scale as an upper limit of the light travel time to the ab- 
sorber. On a direct path, this would place the absorbing medium 
almost 10 times the average orbital diameter of the system away. 
Through scattering in the dense medium, with an optical depth 
T = Nh,i X (Tx ~ 1, the effective light-path is increased, so that 
a size of 2000 It-sec or less is plausible. A sim ilar geometry was 
already proposed bv Mukh eriee & Paull (l2004t) . Through scatter- 
ing in the inner absorber the outer absorber receives a much less 
pulsating X-ray flux and transforms it into fluorescence lines, 
which are consequently also pulsating much less. This effect can 



explain the reduced pulsed fraction in the Fe Ka and Fe KJ3 line 
band (Fig.|7]i. A secondary absorber further out can also explain 
the dip in A^h and covering fraction during the large flare in part 
G (Fig. [TOb . Here the inner absorber is almost completely ion- 
ized and transparent to soft X-rays. The visible absorption is 
then primarily due to the still cold outer absorber. The present 
data do not allow one to distinguish between the lines formed in 
the different absorbers, and as they are both variable it is hard to 
disentangle the effects. The indications we obtain from the data 
do, however, support this absorption geometry. If both columns 
of the partial covering model would be located close to the neu- 
tron star, where the strong column describes the clumpy stellar 
wind and the other the smooth part of the wind, neither the large 
equivalent width nor the strong dip in covering fraction during 
the large flare is easily explainable. The measured widths of the 
Fe Ka and Fe K/3 lines as well as their ratio (Sec. 14.2b can also 
be understood in this geometry as a superposition of cold and 
neutral iron in the outer absorber and higher ionized iron in the 
inner one. As stated before, the spectral resolution of EPIC-pn 
does not allow to separate different iron lines, but we would ex- 
pect to see an only slowly variable Fe Ka line of neutral iron and 
many variable lines from the different ionized species from the 
inner absorber. 



5.2. Fluorescence and abundances 

The 2-10 keV spectrum of GX 301-2 shows evidence for Fe Ka, 
Fe K/3, Ni Ka, Ni K/3, S Ka, CaKa, ArKa, and CrKa fluores- 
cent emission lines. To check consistency with older data, es- 
pecially regarding the CrKa line, we applied our model t o the 
Chandra HETG data presented by (IWatanabe et al.L 12003 1) and 
taken during the pre-periastron flare in January 2001. The model 
describes the data very well in its overall shape. Neither a CrKa 
nor a Ni KJ3 line is significantly detected, but lines with the same 
flux and width as in the XMM-Newton spectram are also not 
ruled out by the data and could vanish due to the worse S/N of 
the Chandra data. 

To describe the XMM-Newton spectrum we used a model in 
which the fluorescent lines were only absorbed by the galactic 
absorption column of A^^^' = 1.7 x 10^^ cm"^ to accommodate 
for the observed fluxes of the Ka lines of S and Ar, see Sect. 14. II 
The flux of the lines indicates, however, that they originate in a 
dense absorber and that thus self-absorption must be taken into 
account. Additionally timing analysis has shown that the main 
producing region of the iron line must be very close to the X-ray 
source and it is likely that this is true for the other lines as well. It 
is thus conceivable that the photons emitted from the fluorescent 
ions are also absorbed by material relatively close to the neutron 
star with a column density in excess of A^^'*', e.g., by the stellar 
wind of the companion. 

In order to estimate how strongly these lines are absorbed 
we performed the following sim ple sanity che ck, assuming so- 
lar abundances for all elements dWilms et all 12000,) : for each 
element, where a fluorescent line was observed, we calculated 
how many photons were absorbed by this element alon e above 
its respective K-shell, applying the cross-sections by Vern er et"al] 
( 1996). Using the fluorescence yields for near neutral ions from 
Kaastra & Mew3 (Il993i) we calculated how many of these pho- 
tons were reemitted in the respective Ka and KJ3 lines. As we 
do not know the geometry of the fluorescent medium and hence 
can not tell how many of these reemitted photons reach us, we 
calculated the fluxes of all lines relative to Fe Ka. We compared 
these relative fluxes to the fluxes as measured by the Gaussian 
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lines in our model. The relative fluxes of the Gaussian lines are 
strongly correlated with the absorption column applied, espe- 
cially for the soft S and Ar lines. By changing the absorption 
column we can find the best agreement between measurement 
and calculation. We find that this is achieved with an overall ab- 
sorption column around 20 x 10^^ cm"^, i.e., a factor of ~6 less 
than the continuum absorption and consistent with the above as- 
sumptions. Systematic uncertainties outweigh any statistical un- 
certainties in this analysis and are estimated to be on the order of 
+5 X 10^^ cm"^. A drastic improvement on this estimate can be 
achieved using Monte Carlo simulations of different absorber ge- 
ometries, but this investigation is beyond the scope of this paper 
(see Barragan et al., 2011, in prep.). 

5.3. The chromium line 

We have reported the detection of an emission line at 5.43 keV, 
compatible with the Ka line energy of neutral chromium. This 
line is weak and only visible because of the large absorption 
column of GX 301-2 making it bright enough to be visible 
over the continuum. The highly absorbed IGR sources described 
above are typically not bright enough to provide enough S/N to 
allow searches for a Cr line. In Supernova Remnants (SNR) evi- 
dence for fluorescent chromium K lines was f ound, most notably 
in W49B dHwang et all 120001: [M iceli et aH l2006h but also in 
Tycho, Cas A, and Kepler dYang et al.i,i200 9x Cr is a decay prod- 
uct of Fe in supernovae, and thus a high iron abundance should 
be accompanied by a high chromium abundance. The abundance 
ratio of Cr to Fe is ~1% for solar mass distribution. In super- 
novae the ratio depends on the explosion model and can be 2% 
or larger for delayed detonation models (see, e.g., Iwamoto et al.', 
Ll999). The SNRs investigated by Yang et al. (2009) show differ- 
ent abundance ratios between 0.63% - 3.6%, leading the authors 
to suspect that they originate from different type of supernovae 
explosions. Using the method described in the previous Section 
to calculate the relative strengths and therewith the relative abun- 
dances of Cr to Fe we obtain a value of 1 .4%, very well within 
the expected values. Taking only the ratio of the measured pho- 
tons in the CiKa and Fe Ka line, we obtain higher values of 
~6.4%, but this value neglects the fact that the K-edge for Cr is 
lower than the one of Fe and that consequently more photons are 
able to ionize Cr With our measurement we can not infer possi- 
ble explosion models or production channels of the GX 301-2 
system. The measured relative line flux does, however, agree 
with predicted fluorescence yields and element abundances. 

5.4. Accretion during the off-state 

The off-state is perhaps the single most puzzling feature of the 
XMM-New ton data. A similar off-state was only recently dis- 
covered by iGo^ii^ et all (l201lh in RXTE data. RXTE's spec- 
tral resolution does not allow for a detailed analysis of the 
iron line region but provi des a larger b r oad ba nd coverage com- 
pared to XMM-Newton. iGogii^ et aP (l201lh fitted individual 
spectra with exposure times of one pulse period each, keeping 
the iron line fixed and letting the photon index of the power 
law continuum vary. They found that the photon index is dras- 
tically increased during the dip, making the overall spectrum 
much softer It thus seems very likely that the accretion was 
switched off during that time and that only weak Comptonization 
in the accretion columns was possible, while most of the X- 
ray flux originated from the black body on the neutron star's 
surface. This result is in accordance with a theory developed 



to explain similar off-state in Vela X-1 (iKrevkenbohm et all 
2008, and references therein). According to this theory, ac- 
cretion ceases during the off-state and the source enters the 
propeller-regime, i.e., a regime i n which the Alvfe n radius is 
larger than the co-rotating radius dLamb et al.Lll973l) . This tran- 
sition is triggered by a strong decrease in the mass accretion 
rate, possibly d ue to a strongly structured or clumpe d wind 
(Oskinov ae t al l 120071: Krevkenbohmet al., 2008: Fur st et all 
201Q). Even though GX 301-2 is not directly accreting from 
the wind in th e pre-p eriastron flare b ut rather from an ac- 
cretion stream dLeahv & Kostkai |2008|) . strong density varia- 
tions are evident through the flaring behavior of the light-curve 
(Fig.|2]l. Other observational evi dence for a clumpy wind a round 
Wray 977 was put forward by iMukheriee & Paull d2004t) who 
studied in detail the time variability of the absorption column on 
short and long timescales. 

All our results indicate that the observed off-state was caused 
by a cessation of accretion and not by an optic ally- thick clump 
moving through the line of sight. A clump in the line of sight 
would also lead to drastically reduced pulsations and flux, as 
the X-ray flux is scattered in the absorber and thereby smearing 
out the pulsations dMakino et al.l Il985h . In the spectral analy- 
sis, however, we would expect that the spectrum gets harder, as 
softer X-rays are affected stronger by absorption than harder X- 
rays. We would also expect to see a drastic increase of the ab- 
sorption column. Both are not the case, as seen in Fig.fTOl 

We investigated the overall behavior of the pulsed fraction as 
function of A^h, and found almos t no correlati on (Fig. [13]), con- 
tradictory to the results by Makin oet al.l dl985h . Scattering in the 
absorption medium seems to be influencing the pulsed flux only 
marginally, so that the complete disappearance of pulses in part 
D is very unlikely to be caused by larger A^h values. The variation 
in the pulsed fraction seems to originate from different accretion 
rates, as expected from a clumpy accretion medium. This influ- 
ence is also supported by the variability of the individual pulse 
profiles, as shown in Fig. [3] We do not see a weaker pulsed flux 
with essentially the same shape causing the lower pulsed frac- 
tion, but rather a clear change in the overall pulse profile caused 
by a different accretion rate. 

From results at harder X-rays dEvangelista et al.Ll20ld) a de- 
crease of the pulsed fraction with decreased luminosity would 
be expected. In a simple toy model this can be explained by 
a change in altitude of the X-ray producing region in the ac- 
cretion column (Lutovinov & Tsygankov, 2009). Such a change 
would also be visible in a change of the en ergy of the CRSF 
dSuchv et all 1201 it iDoroshenko et al.L l2010l) . We did not find 
a correlation between the pulsed flux and the luminosity in 
the XMM-Newton data. The pulsed fraction during the pre- 
periastron flare is the lowest over the whole orbit and is also 
much lower in the soft than in the hard X-rays dEvangelista et all 
2010). It is thus possible that a change in the accretion column 
might not be visible in the pulsed fraction of the 2-lOkeV band. 

Assuming that the off-state was caused by a cessation of ac- 
cretion, possibly because the neutron star entered the propeller 
regime, the remaining X-ray flux is mainly coming from the 
still fluorescent absorbing medium around the neutron star, as 
described in Sec. 15. II 

As shown in Fig. [T2] during the normal states the equivalent 
width of the Fe Ka is correlated with the A^h.i value. This correla- 
tion is explained by the assumption that the fluorescent medium 
is the same as the absorbing medium. As a result, a higher ab- 
sorption column leads also to a brighter fluorescence line as 
more atoms in the absorber are ionized by the X-ray continuum. 
During the off-state, the X-ray continuum is drastically reduced 



12 



F. Flirst et al.: Fluorescent lines and absorption in GX 301-2 



1 

0.9 
0.8 
0.7 



§ 0.1 



I 0.5 

T3 

0) 

0.4 

a 



0.3 
0.2 
0.1 



- 




100 120 140 160 180 
A^H.i [10^2 H-atoms cm- 



200 



220 



Fig. 13. Pulsed fraction as function of A^h.i- Error bars on the 
pulsed fraction are arbitrary. Color code is the same as in Fig.fTOl 



leading to a reduced Fe Ka line from the main absorber in the 
line of sight. The off-state indicates a very clumpy medium, so 
that the absorption in the line of sight is not necessarily an in- 
dicator for the overall density of the medium surrounding the 
X-ray source. The observed break-down of the correlation be- 
tween equivalent width of the Fe Ka line and the A^h.i value can 
thus be explained by a complex structure of the medium around 
the neutron star For example, fluorescence from dense medium 
behind the neutron star could arrive simultaneously with the re- 
duced continuum flux at the observe due to different light-travel 
times. The remaining continuum flux during the off-state could 
also be due to Compton reflection from a dense mediu m around 
the n eutron star, similar to the scenario proposed by iRea et al] 
(1200 5) to explain a long low-intensity state in the Low Mass 
X-ray Binary (LMXB) GX 4+1. These authors analysed broad 
band BeppoSAX spectra and found that they could be described 
by a cut-off power law plus a reflection component, the latter 
becoming more important during the low intensity state. In this 
state also the Fe Ka l ine reached the hightes equivalent widths, 
similar to GX 301-2. iRea et al.1 (12005.) found that although the 
pulsations vanished below 7 keV they could still be seen in the 
harder X-rays, an effect we can not investigate with the XMM- 
Newton data but wh ich is not evident in the RXTE/PCA data 
presented bv iGogii^ et aP (l201 Ih . If the reflecting surface is not 
as well defined as in the geometry proposed for GX 1 H-4, all pul- 
sations can be smeared out over a broad energy range, so that 
the reflected continuum flux is also unpulsed when escaping the 
medium around the system. More observations in a broader en- 
ergy range are necessary to understand the geometry further. 

5.5. Outlook 

Off-state are seen in an increasingly number of sources, like 
Vela X-1 and GX 1+4 and can provide interesting insight in the 
structure of the accreted material as well as the configuration 
of the accretion and absorption region. For GX 301-2 we pro- 
pose that two absorbers exist, a very dense but clumpy primary 
one close to the neutron star, and a weaker and smoother sec- 



ondary one farther away. Fluorescence lines from the secondary 
absorber are much weaker then the ones from the primary and 
only visible during the off-state. High-resolution spectra of the 
off-states from Chandra HETG or future missions will help to 
investigate them in more detail by allowing for better measure- 
ments of the ionization state and the geometry of the absorber 
The discovery of a Cr Ka line in the spectrum of GX 301-2 
leads to interesting new questions about the abundances of the 
iron family atoms, including Cr and Ni. With higher effective 
areas in future missions, more sources might be found show- 
ing a Cr Ka line and in comparing these findings to GX 301-2 
more insight into the production of HMXB and the preceding 
supernova can be achieved. With broad-band spectroscopy, as al- 
ready possible with the Suzaku observatory, the underlying con- 
tinuum can be constrained simultaneously in a better way than 
with XMM-Newton oder Chandra only. Our collaboration is an- 
alyzing Suzaku data of GX 301 -2 taken dur ing orbital phases 
outside the pre-periastron flare (iSuchv et al.L [201 1) . These ob- 
servations provide insight into the behavior of the accretion col- 
umn through exact measurements of the variability of the CRSF 
energy. All this information combined will complete the picture 
of GX 301-2 and possibly also of the other highly absorbed 
sources. 
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